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Abstract. Radial pulsations of newborn neutron stars 
(protoneutron stars) are studied for a range of internal 
temperatures and entropies per baryon predicted by the 
existing numerical simulations. Protoneutron star models 
are constructed using a realistic equation of state of hot 
dense matter, and under various assumptions concerning 
stellar interior (large trapped lepton number, zero trapped 
lepton number, isentropic, isothermal). Under prevailing 
conditions, linear oscillations of a protoneutron star can 
be treated as adiabatic, and evolutionary effects can be 
neglected on dynamic timescale. The dynamic behavior is 
governed by the adiabatic index, which in turn depends 
on the physical state of the stellar interior. The eigenfre- 
quencies of the lowest radial modes of linear, adiabatic 
pulsations are calculated. Stability of protoneutron stars 
with respect to small radial perturbations is studied, and 
the validity of the static stability criteria is discussed. 
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1. Introduction 

Newly born neutron stars are hot and lepton rich objects, 
quite different from ordinary low temperature, lepton poor 
neutron stars. In view of these differences, newly born neu- 
tron stars are called protoneutron stars; they transform 
into standard neutron stars on a timescale of the order 
of ten seconds, needed for the loss of a significant lepton 
number excess via emission of neutrinos trapped in dense, 
hot interior. 

In view of the fact that the typical evolution timescale 
of a protoneutron star (seconds) is some three orders of 
magnitude longer, than the dynamical timescale for this 
objects (milliseconds), one can study its evolution in the 
quasistatic approximation (Burrows & Lattimer 1986). 
Static properties of protoneutron stars, under various as- 
sumptions concerning composition and equation of state 
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(EOS) of hot, dense stellar interior were studied by numer- 
ous authors (Burrows & Lattimer 1986, Takatsuka 1995, 
Bombaci et al. 1995, Bombaci 1996, Prakash et al. 1997). 

The scenario of transformation of a protoneutron star 
into a neutron star could be strongly influenced by a phase 
transition in the central region of the star. Brown and 
Bethe (1994) suggested a phase transition implied by the 
K'^ condensation at supranuclear densities. Such a 
condensation could dramatically soften the equation of 
state of dense matter, leading to a low maximum allowable 
mass of neutron stars. In such a case, the massive pro- 
toneutron stars could be stabilized by the effects of high 
temperature and of the presence of trapped neutrinos, and 
this would lead to maximum baryon mass of protoneutron 
stars larger by some 0.2 Mq than that of cold neutron 
stars. The deleptonization and cooling of protoneutron 
stars of baryon mass exceeding the maximum allowable 
baryon mass for neutron stars, would then inevitably lead 
to their collapse into black holes. The dynamics of such a 
process was recently studied by Baumgarte et al. (1996). It 
should be mentioned, however, that the very possibility of 
existence of the kaon condensate (or other exotic phases of 
matter, such as the pion condensate, or the quark matter) 
at neutron star densities is far from being established. Re- 
cently, for instance, Pandharipande et al. (1995) pointed 
out, that kaon-nucleon and nucleon-nucleon correlations 
in dense matter raise significantly the threshold density 
for kaon condensation. In view of these uncertainties, we 
will restrict in the present paper to a standard model of 
dense matter, composed of nucleons and leptons. 

The calculations of the static models of protoneutron 
stars should be considered as a first step in the studies of 
these objects. It is clear, in view of the dynamical scenario 
of their formation, that protoneutron stars are far from 
being static. The formation scenario involves compression 
(with overshoot of central density) and a hydrodynamical 
bounce, so that a newborn protoneutron star begins its life 
in a highly excited state, pulsating around its quasistatic 
equilibrium. If the rotation of protoneutron star is not 
too rapid, the coupling of radial and non-radial pulsations 
is weak. Such a scenario of formation leads to a prefer- 
ential excitation of the radial modes. Having constructed 
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the static model of a protoneutron star, one should thus 

answer the qiiestion about the stability of the static con- 
figuration with respect to radial perturbations (standard 
stability criteria are valid only for cold neutron stars). 
Clearly, both high temperature and large trapped lepton 
number will influence the spectrum of radial eigenfrequen- 
cies of protoneutron stars, implying differences with re- 
spect to the case of cold neutron stars. 

In the present paper we study the radial pulsations of 
protoneutron stars and their stability. Our models of pro- 
toneutron stars are composed of a hot, neutrino-opaque 
interior (hereafter referred to as "hot interior"), sepa- 
rated from much colder, neutrino-transparent envelope 
by a neutrinosphere. We will consider two limiting cases 
of the thermal state of the hot interior: isentropic, with 
entropy per baryon s = const., and isothermal, with 
Tqc = [gaaY^'^T = const. [Too is the value of the temper- 
ature, measured by an observer at infinity, while T is the 
value of the temperature measured by a local observer). 

The first case, characteristic of a very initial state of 
a protoneutron star, will simultaneously correspond to a 
significant trapped lepton number. The second case cor- 
responds to situation after the deleptonization of a pro- 
toneutron star. The position of the neutrinosphere will be 
located using a simple prescription based on specific prop- 
erties of the neutrino opacity of hot dense matter. In all 
cases, the equation of state of hot dense matter will be 
determined using one of the realistic models of Lattimer 
and Swesty (1991). 

The plan of the paper is as follows. In Section 2 we 
describe the physical state of the interior of protoneutron 
star, with particular emphasis on the EOS of the hot inte- 
rior at various stages of evolution of a protoneutron star. 
We explain also our prescription for locating the neutri- 
nosphere of a protoneutron star, and we give some de- 
tails concerning the assumed temperature profile within 
protoneutron star. Equations of state, corresponding to 
different physical situations, are described in Section 3. 
We also present there static models of protoneutron stars, 
calculated for various assumptions concerning the hot stel- 
lar interior. In Section 4 we compare various timescales, 
characteristic of evolution and dynamics of a protoneu- 
tron star, which are essential for approximations used in 
the treatment of pulsations. Section 5 is devoted to the for- 
mulation of the problem of linear, adiabatic, radial pulsa- 
tions of protoneutron stars. Both pulsations and stability 
involve the adiabatic index of pulsating stellar interior; 
our results for this important quantity are presented in 
Section 6. Numerical results for the eigenfrequencies of 
the lowest modes of radial pulsations and the problem of 
stability of protoneutron stars are presented in Section 7. 
Finally, Section 8 contains a discussion of our results and 
conclusion. 



2. Physical model of the interior of protoneutron 
stars 

We will consider a protoneutron star (PNS) just after its 
formation. We will assume it has a well defined neutri- 
nosphere, which separates a hot, neutrino-opaque interior 
from colder, neutrino-transparent outer envelope. The pa- 
rameters, which determine the local state of the matter 
in the hot interior are: baryon (nucleon) number density 
n, net electron fraction Y,, = («.,.- — n,.j-)/n, and the net 
electron- neutrino fraction — Y^^ — Yp^. The calcula- 
tion of the composition of hot matter and of its EOS is 
described below. 

2.1. Neutrino opaque core with nonzero trapped lepton 
number 

The situation described in this subsection is characteristic 
of the very initial stage of existence of a PNS. Matter is as- 
sumed to be composed of nucleons (both free and bound in 
nuclei) and leptons (electrons and neutrinos; for the sake 
of simplicity, we do not include muons). All constituents 
of the matter (plus photons) are in thermodynamic equi- 
librium at given values of n, T and Yi = Y^ + Y,^. The 
composition of the matter is calculated from the condi- 
tion of beta equilibrium, combined with the condition of 
a fixed Yi , 

Yi = Y, + Y,, (1) 

where fij are the chemical potentials of matter con- 
stituents. At the very initial stage we expect Yi ~ 0.4. 
Electron neutrinos are degenerate, with fi^^ ^ T (in what 
follows we measure T in energy units). The deleptoniza- 
tion, implying the decrease of V; , occurs due to diffusion 
of neutrinos outward (driven by the fi^^ gradient), on a 
timescale of seconds (Sawyer & Soni 1979, Bombaci et al. 
1996). The diffusion of highly degenerate neutrinos from 
the central c;ore is a dissipative process, resulting in a sig- 
nificant heating of the neutrino-opaque core (Burrows & 
Lattimer 1986). 

2.2. Neutrino opaque core with Y^ =0 

This is the limiting case, reached after complete delep- 
tonization. There is no trapped lepton number, so that 
Yi = Yf, and Y^^ = Yp^, and therefore fj.„^ = Hi?e = 0. 
Neutrinos trapped within the hot interior do not there- 
fore influence the beta equilibrium of nucleons, electrons 
and positrons, and for given n and T the equilibrium value 
of Ye is determined from 

l^p + lJ-e = fJ-n , (2) 

while /ie+ = — /ig. In practice, this approximation can be 
used as soon as electron neutrinos become non-degenerate 
within the opaque core, /i^^ < T, which occurs after some 
~ 10 seconds (Sawyer & Soni 1979, Prakash et al. 1997). 
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2.3. Neutrino sphere and the temperature profile 

In principle, the temperature (or entropy per nucleon) 
profile within a PNS has to be determined via evolution- 
ary calculation, starting from some initial state, and tak- 
ing into account relevant transport processes in the PNS 
interior, as well as neutrino emission from PNS. Trans- 
port processes within neutrino-opaque interior occur on 
timescales of seconds, some three orders of magnitude 
longer than dynamical timescales. Convection can shorten 
neutrino transport timescale, but still the time needed for 
the deleptonization of the neutrino opaque core is then 
much longer than the dynamical timescale. The very outer 
layer of a PNS becomes rapidly transparent to neutrinos, 
deleptonizes, and cools on a very short timescale via e~e^ 
pair annihilation and plasmon decay to T < 1 MeV. It 
seemed thus natural to model the thermal structure of the 
PNS interior by a hot core limited by a neutrinosphere, 
and a much cooler, neutrino transparent outer envelope. 
The transition through the neutrinosphere is accompanied 
by a temperature drop, which takes place over some inter- 
val of density just above the "edge" of the hot neutrino- 
opaque core, situated at some n^. 

In view of the uncertainties in the actual tempera- 
ture profiles within the hot interior of PNS, we consid- 
ered two extremal situations for n > rin, correspond- 
ing to an isentropic and an isothermal hot interior. In 
the first case, hot interior was characterized by a con- 
stant entropy per baryon s — const.. In the case of 
trapped lepton number, this leads to the EOS of the type: 
pressure P — P{n, [s,Yi]), energy density divided by 
(energy-mass density), p — pin, and temperature 

T = r(n, [s,Yi\), with fixed s and Y;. This EOS will be 
denoted by EOS[s, F/]. 

The condition of isothermahty, which in the static case 
corresponds to a vanishing heat flux, is more complicated. 
Due to the curvature of the space-time within PNS, the 
condition of isothermahty (thermal equilibration) corre- 
sponds to the constancy of T^o = {gooY^'^T (see, e.g., 
Zeldovich & Novikov 1971, chapter 10.6). Actually, the 
isothermal state within the hot interior will be reached on 
a timescale corresponding to thermal equilibration, which 
is much longer than the lifetime of a PNS. Nevertheless, 
we considered the Too = const, models because, as a lim- 
iting case so different from the s = const, one, it enables 
us to check the dependence of our results for PNS on the 
thermal state of the hot interior. 

To determine the isothermal temperature profile in the 
hot interior after deleptonization, we use the condition of 
thermal equilibrium, given by the constancy of the func- 
tion Too = T{r)e^'^^^^ (where ^{r) is the metric fuction, 
see Sect. 5). The relativistic condition of the isothermahty 
can be rewritten as: 

dlnT _ 1 dP 
dr " pc2 + P d^ ■ ^ ^ 



This formula enables us to determine the T(n) profile, for 
given EOS, independently of the specific structure of the 
stellar model under consideration. 

Treating n, T as thermodynamic variables for our 
equation of state, we can rewrite Eq. (^ in the form: 

dlnT _ P 
d In n pc^ + P 

fd]nP\ r _ P f dinP \\-' 

^ \dlnn)\ pc^ + P\d\nT)} ' 

Using the above formula we can construct a specific 
isothermal EOS, describing hot isothermal interiors, and 
parametrized by the boundary condition at the edge of the 
hot isothermal core - the value of T just below the neu- 
trinosphere will be denoted by Tb. Thus in our relativistic 
calculations the set of the "isothermal" stellar configura- 
tions corresponds to stars with given Tb. Starting from 
T — Th, the temperature increases inward, reaching its 
maximum value in the center of the star where T = Tccntr- 
This central temperature is the function of the central den- 
sity pcontr and is larger for a star with larger Pcontr- 

Our calculation of the neutrinosphere within the hot 
PNS interior is explained below. For a given static PNS 
model, the neutrinosphere radius, Ru, has been located 
through the condition 

where is calculated at the matter temperature. En is 
the mean energy of non-degenerate neutrinos at and above 
the neutrinosphere, E^, = 3.15T^, and rpmp is the proper 
distance from the star center. We assumed that neutrino 
opacity above is dominated by the elastic scattering off 
nuclei and nucleons, so that Aj, = X°{n,T)/El. Then, we 
determined the value of the density at the neutrinosphere, 
n^, for a given static PNS model, combining Eq. (||) with 
that of hydrostatic equilibrium, and readjusting accord- 
ingly the temperature profile within the outer layers of 
PNS. 

Neutrino opacity in the outer layers of PNS can be 
well approximated by l/Xi, ~ kqEJ^ . Within a reasonable 
approximation kq ~ 7p, where p is the matter density, 
and 7 = 6- 10"^" cm~^ (p in g cm""^ and E^ in MeV). 
The proper distance near the neutrinosphere is given by 
drprop = [grrY'^dr ~ (1 - 2GM/Pc2)-i/2dr. Using the 
definition of the neutrinosphere radius, R^^ we can express 
the mass of the envelope above R^, and the pressure at R^, 
denoted by P^, in terms of E^ = 3.15T^. This leads to 

^ / m\ ^ era; 

Pn = 2.23 10^2 MR-^ 1 - 0.295^ T^"^ . (6) 
\ R J cm3 

where M = M/Mq and R = R/10 km. 
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For a given stellar model, this approximate relation, 
combined with Eq.(3), enables one to determine, in a 
self-consistent way, the values of T, n, and p at the neu- 
trinosphere. In practice, this was done assuming a spe- 
cific functional form of the temperature drop within the 
neutrinosphere (a combination of Fermi functions), which 
yielded a smooth transition between hot interior and cool 
envelope. In all cases, the temperature profile was adjusted 
in such a way, that neutrinosphere was found within the 
region of the temperature drop. For a 1.4 Mq PNS, and 
an isothermal hot core with Tb = 15 McV, we found Tj, = 
4.3 MeV, n„ = 2.0 lO'^ fm'^, and = 3.5 lO^^ g ^^-3_ 
In the case of a very massive PNS with M = 2 Mq 
we obtained (for the same value of Tb) Ti, = 5.2 MeV, 

= 2.4 10"^ fm"3, and p^ = 4.1 10^^ g cm'^. These re- 
sults are in a reasonable agreement with those obtained in 
detailed numerical simulations (see, e.g.. Burrows & Lat- 
timer 1986). 

3. Equation of state and static models of protoneu- 
tron stars 

The starting point for the construction of our EOS for 
the PNS models was the model of hot dense matter of 
Lattimer and Swesty (1991), hereafter referred to as LS. 
Actually, we used one specific LS model, corresponding to 
the incompressibility modulus at the saturation density of 
symmetric nuclear matter K = 220 MeV (this model will 
be hereafter referred to as LS-220). For n > n,y we sup- 
plemented the LS-220 model with contributions resulting 
from the presence of trapped neutrinos of three flavours 
(electronic, muonic and tauonic) and of the corresponding 
antineutrinos. 

In Figs. 1-3 we show our EOS in several cases, cor- 
responding to various physical condition in the hot, 
neutrino-opaque interior of PNS. For the sake of com- 
parison, we have also shown the EOS for cold catalyzed 
matter, used for the calculation of the (cold) NS models. 
In Fig. 1 we show our EOS at subnuclear densities. At 
such densities, both the temperature and the presence of 
trapped neutrinos lead to a significant increase of pressure, 
as compared to the cold catalyzed matter. The constant 
Too EOS stiffens considerably at lower densities, which is 
due to the weak dependence of the thermal contribution 
(photons, neutrinos) on the baryon density of the mat- 
ter. It is quite obvious, that Too = const. EOS becomes 
dominated by thermal effects for n < 10~^ fm~^. On the 
contrary, for the isentropic EOS, the effect of the trapped 
lepton number (Yj = 0.4) turns out to be more important 
than the thermal effects. This can be seen in Fig. 1, by 
comparing long-dashed curve, [s = 2, Yi = 0.4], with a 
short-dashed line, which corresponds to an artificial (un- 
physical) case with small thermal effects, [s = 0.5, Yi = 
0.4]. 

It is clear, that the correct location of the neutri- 
nosphere, which separates hot interior from the colder 
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Fig. 1. Pressure versus baryon density for our model of dense 
hot matter (based on the LS-220 model of the nucleon com- 
ponent of the matter), under various physical conditions, for 
the subnuclear densities n < no = 0.16 fm^"^. The curve T — 
corresponds to cold catalyzed matter. The curve correspond- 
ing to s = 0.5, Yi = 0.4 is unphysical, but has been added in 
order to visualize the importance of trapped lepton number at 
subnuclear densities. The low density edge of the hot, neutrino 
opaque core corresponds to n-i^ = 2 x 10~^ fm~^. 



outer envelope, should be important for the determina- 
tion of the radius of PNS. 

Our EOS above nuclear density is plotted in Fig. 2. 
The presence of a trapped lepton number softens the EOS, 
while thermal effects always lead to pressure increase with 
respect to that for cold catalyzed matter. The softening of 
the supranuclear EOS at fixed Yi is due to the fact, that 
a significant trapped lepton number increases the proton 
fraction, which implies the softening of the nucleon con- 
tribution to the EOS. This is the reason why pressure for 
[s — 0.5, Yi = 0.4] model is smaller than in the case of 
cold, catalyzed matter. 

In the calculations of the stellar structure the rele- 
vant EOS is of the form P = P{p), because only pressure 
and mass-energy density appear in the general relativis- 
tic equations of hydrostatic equilibrium. It is obvious that 
thermal contribution to p is always positive. Also, the con- 
tribution of trapped neutrinos (which do not contribute 
to baryon number density n) to p is positive. However, 
large trapped lepton number implies an increased proton 
fraction, which in turn implies a softening of the nucleon 
contribution to the pressure. The interplay of these effects 
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Fig. 2. Pressure versus baryon density for our model of dense 
hot matter, under various physical conditions, for the supranu- 
clear densities. The curve T = corresponds to cold catalyzed 
matter. The curve corresponding to s = 0.5, Yi = 0.4 is un- 
physical, but has been added in order to visualize the impor- 
tance of trapped lepton number at subnuclear densities. The 
low density edge of the hot, neutrino opaque core corresponds 
ton^ = 2x 10"^ fm"^. 



leads to a characteristic softening of the P(p) EOS for 
large trapped lepton number, visualized in Fig. 3. 

It should be stressed, that in contrast to Bombaci et 
al. (1995) and Prakash et al. (1997) we used a unified 
dense matter model, valid for both supranuclear and sub- 
nuclear densities. Also, the fact that we use various as- 
sumptions about the T and s profiles within PNS, enables 
us to study the relative importance of the temperature 
profile and that of a trapped lepton number, for the PNS 
models. 

The mass-radius relation for the PNS models calcu- 
lated using various versions of our EOS for the hot inte- 
rior is shown in Fig. 4. We assumed = 2 x 10~^ fm~^, 
which was consistent with our definition of the ncutri- 
nosphere. For the sake of comparison, we also show the 
mass-radius relation for the T = (cold catalyzed mat- 
ter) EOS, which corresponds to cold neutron star models. 
In the case of an isothermal hot interior with T\y = 15 MeV 
we note a very small increase of the maximum mass, as 
compared to the T = case (this result is consistent with 
those obtained by Bombaci et al. 1995, Prakash et al. 1997, 
for models composed of nucleons and leptons only). How- 
ever, the effect on the mass-radius relation is quite strong. 



Fig. 3. Pressure versus mass-energy density for our model of 
dense hot matter, under various physical conditions. The curve 
T = corresponds to cold catalyzed matter. 

and increases rapidly with decreasing stellar mass. In the 
case of the isentropic EOS with a trapped lepton number, 
[s ~ 2,Yi ~ 0.4], the softening of the high-density EOS 
due to the trapped Yi leads to the decrease of Mmax com- 
pared to the T = far as the value of Mmax is 
concerned, the softening ciT{x;t of Yi prevails over that of 
finite s (this is consistent with results of Takatsuka 1995, 
Bombaci et al. 1995, and Prakash et al. 1997 for purely 
micleonic models). However, the thermal effect on the stel- 
lar radius is important even in the case of 1; = 0.4. This 
can be seen by comparing [s = 2, Yi = 0.4] curve with 
that corresponding to the unphysical, fictitious case of 
[s = 0.5, Yi = 0.4]. Let us notice, that for our EOS, con- 
figurations with maximum (minimum) gravitational mass 
are simultaneously those with the maximum (minimum) 
baryon mass. The values of the minimum and maximum 
gravitational and baryon masses for various EOS of the 
PNS interiors, based on the LS-220 models of the mi- 
cleon component of dense matter, are given in Table 1. 
The baryon (rest) mass of a PNS model is defined by 
Mbiir = Am,o, where A is the total baryon number of a 
PNS, and mo is the mass of the hydrogen atom. 

We find that the value of Mbar.max is the largest one 
for the T = (cold catalyzed matter) EOS. For an isother- 
mal hot star with Tb = 15 MeV the total baryon mass is 
slightly lower than M^~^^^, but for isentropic stars we 
see a significant decrease of Mbar.max as compared to the 
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Fig. 4. The gravitational mass versus stellar radius for static 
models of the protoneutron stars and neutron stars, under var- 
ious assumptions concerning the physical conditions within the 
stellar interior. The curve corresponding to s = 0.5, Yj = 0.4 
is unphysical, but has been added in order to visualize the rel- 
ative importance of the trapped lepton number and thermal 
effects. The curve T = corresponds to cold neutron stars. 



Fig. 5. The gravitational mass versus central density for static 
models of the protoneutron stars and neutron stars, under var- 
ious assumptions concerning the physical conditions within the 
stellar interior. The curve corresponding to s = 0.5, Yi = 0.4 
is unphysical, but has been added in order to visualize the rel- 
ative importance of the trapped lepton number and thermal 
effects. The curve T = corresponds to cold neutron stars. 



Table 1. Minimum and maximum gravitational and baryon 
masses of cold neutron star (T = 0) and isothermal 
(Tb = 15 MeV, = 0) and isentropic (s = 2,Yi = 0.4) pro- 
toneutron stars. Equations of state are based on the LS-220 
model of the nucleon component of the matter. 



EOS 



Mk 



[Mq] [Me 



[Me] [Me] 



T = 0, = 0.054 0.055 2.044 2.406 

T = 15, n = 0.864 0.892 2.050 2.391 
s = 2,Yi= 0.4 0.675 0.676 1.972 2.183 



T = case. A newborn hot protoneutron star with the 

maximum gravitational mass mL=x^^^-°-^1 = 1.97 Me 
transforms, due to deleptonization and cooling, into a cold 
NS of gravitational mass of 1.9 Mq, some 0.15 Mq less 
than the maximum gravitational mass of cold neutron 
stars (c.f. the analysis of Bombaci 1996). 

For a given model of the neutrino opaque PNS core, 
the value of Mmin depends on the location of the neutri- 
nosphere. However, this dependence is relatively weak. For 



example in the case of EOS [s = 2, Yi — 0.4] a rather dras- 
tic variation of rii, from 2 x 10~^ fm""^ to 6 x 10^"^ fm"*^ im- 
plies change in mK^^' from 0.686 Mq to 0.649 Mq, 
respectively (the radius of these minimum-mass stars 
varies then from 44 km to 32.8 km, respectively). 

At a given mass, the radius of a PNS is significantly 
larger than that of a cold NS. It should be stressed, how- 
ever, that the value of radius turns out to be quite sensitive 
to the location of the edge of the hot neutrino-opaque in- 
terior (i.e., to the value of n^), especially for PNS which 
are not close to the Mmax configuration. The choice of 
Prakash et al. (1997) would lead to a much smaller effect 
on R. 

The role of the thermal and the lepton number effects 
are particularly pronounced in the mass - central density 
plots for the PNS models, Fig. 5. The low - Pcentr segments 
of the M — pcontr plots for the PNS arc dramatically dif- 
ferent from that for the cold NS models. The relevance of 
this features for the stability of lower-mass PNS will be 
discussed in Sect. 7. 
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4. Characteristic timescales 

The EOS of PNS is evolving with time, due mainly to 
the deleptonization process, which changes the composi- 
tion of the matter, and also due to changes of the inter- 
nal temperature of the star. However, these changes oc- 
cur on the timescales Tevoi ~1-10 s (see below), which are 
three or more orders of magnitude longer than the dy- 
namical timcscale, governing the readjustment of pressure 
and gravity forces. This dynamical timescale r^yn ^ 1 nis 
corresponds also to the characteristic periods of the PNS 
pulsations. In view of this, we are able to decouple the PNS 
evolution from its dynamics, with a well defined EOS of 
the PNS matter. 

The evolution of the PNS interior results from cool- 
ing, deleptonization, and dissipative transport processes, 
leading to heating. Deleptonization of PNS is due to dif- 
fusion of Uf,, driven by the gradient of their chemical po- 
tential, and occurs on a timescale rjeiept ~ few seconds 
(Prakash et al. 1997). (Both Tdeicpt and TcooI can actually 
be shorter, due to the presence of convection in some lay- 
ers of a PNS). The PNS cooling is due to neutrino losses 
from the neutrinospherc; for the neutrino-opaque interior, 
the characteristic timescale of cooling TcqoI is of the or- 
der of tens of seconds (Sawyer & Soni 1979). Notice, that 
deleptonization is accompanied by a significant heating 
of neutrino-opaque core. Highly degenerate u^, diffuse out 
from the neutrino-opaque core, and due to E,^^ » T they 
deposit most of their energy in the matter, which in view 
of Tcooi ^ Tdeiept Corresponds to the net heating. 

Radial pulsations of PNS are damped due to dis- 
sipative processes, resulting from the weak interactions 
involving nucleons and leptons. The dissipative effects 
can be represented in the form of a bulk viscosity of 
hot dense matter (Sawyer 1980). Let us consider first 
the case of hot interior with a significant trapped lep- 
ton number. For Yi = 0.3 Sawyer (1980) gets Tdamp{Yi = 
0.3) ~ 2000 (T/10 MeV)2 s. For deleptonized hot in- 
terior the characteristic timescale is somewhat shorter, 
Tdainp(>^ = 0) ~ 30 (T/10 MeV)2 s (Sawyer 1980). Still, 
in both cases we get Tdamp>10"'Tdyn, so that one can safely 
neglect damping when calculating the eigenfrequencies of 
radial pulsations of PNS. 

Summarizing, the estimates of the evolutionary and 
dissipative timescales indicate, that radial pulsations of 
the hot interior of PNS can be treated as adiabatic, and 
can be studied using a well defined EOS of dense hot mat- 
ter, corresponding to a given stage of evolution of a PNS. 

Of course, all these remarks are valid only for the hot 
interior of PNS (i.e., the region below the neutrinosphere). 
However, the outer envelope contains less than 10"^ of the 
mass of PNS, and its infiuence on the eigenfrequencies of 
radial pulsations is negligible. 



5. Linear adiabatic radial pulsations of protoneu- 
tron stars 

Consider an idealized static configuration of a PNS and 
assume it is spherically symmetric. Using the notation of 
Landau & Lifshitz (1975), we write the metric for such a 
configuration as 



ds' = e^c'dt' - eMr' - r'ide' + sin^ M</)^) , 



(7) 



where A and are functions of r. 

The hydrostatic equilibrium of the static PNS is de- 
scribed by the Tolman-Oppenheimer-Volkoff (TOV) equa- 
tions (Tolman 1939, Oppenheimcr & Volkoff 1939) 



dP 
dr 

dm 
dr 

diy 
dr 



Gmp 



A-Kr^p 



dP 



1 



P 



1 + ^1 , (8) 



(9) 



(P + pc^) dr 



(10) 



where m is the mass contained within radius r. Since our 
EOS of PNS can be always written in the one-parameter 
form P = P{p), the TOV equations can be numerically 
integrated for a given central density Pcentrj yielding the 
stellar radius, R, and the total gravitational mass, M = 
m{R), of the star. 

The equations governing infinitesimal radial adiabatic 
stellar pulsations in general relativity were derived by 
Chandrasekhar (1964), and were rewritten by Chanmu- 
gan (1977) in a form, which turns out to be particularly 
suitable for numerical applications. Two important quan- 
tities, describing pulsations, are: the relative radial dis- 
placement, ^ = Ar/r, where Ar is the radial displace- 
ment of a matter element, and AP - the corresponding 
Lagrangian perturbation of the pressure. These two quan- 
tities are determined from a system of two ordinary dif- 
ferential equations, which we rewrite as 



^ = 1 f^, AP\ dP e 
dr r V TP y dr {P + pc^) 



(11) 



dAP X 2\ AP 



dr 



'dP 



87rG 



dr J {P + pc^) c4 



e^{P + pc^)Pr\ 



+ AP 



A-kG . 



dP 1 



dr (P + pc^) 



(P + pc2)reM ,(12) 



where P is a relativistic adiabatic index (see Section 6), 
LO is the eigenfrequency and the quantities ^ and AP are 
assumed to have a harmonic time dependence oc e*'^*. Our 
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Eq. (12) has been obtained from a second-order pulsa- 
tion equation of Chandrasekhar (1964) [his Eq. (59)], us- 
ing his Eq. (35) and Eq.(36). While Eq. (11) coincides with 
Eq.(18) of Chanmugan (1977), our Eq. (12) - in contrast 
to his second equation [Eq.(19) of Chanmugan 1977] - does 
not show a singularity at the stellar surface. Another im- 
portant advantage of our system of pulsation equations, 
Eq. (11, 12), stems from the fact, that they do not involve 
any derivatives of the adiabatic index, T. In view of the 
fact, that we used tabulated forms of the EOS and of F, 
this enabled us to reach a very high precision in solving 
the eigenvalue problem. 

To solve equations (11) and (12) one needs two bound- 
ary conditions. The condition of regularity at r = re- 
quires, that for r ^ the coefficient of the 1/r-term in 
Eq. (11) must vanish. 



(AP), 



-3(^rp)^ 



(13) 



Our normalization of eigenfunctions corresponds to ^(0) — 
1. The surface of the star is determined by the condition 
that for r ^ R, one has P ^ 0. This implies 



(AP) 



surface 



= 



(14) 



The problem of solving Eqs. (11), (12), can be reduced 
to a second order (in ^) , linear radial wave equation, of the 
Sturni-Liouville type. The quantity is the eigenvalue 
of the Sturm-Liouville problem with boundary conditions 
given by Eqs. ([T^, (p^). For a given static model of a PNS, 
we get a set of the eigenvalues iv^ < ujf < ... < uj^ < 
with corresponding eigenfunctions ^o, ^i, ^n, where 
the eigenfunction has n nodes within the star, < r < 
R (see, e.g.. Cox 1980). 

A given static model is stable with respect to small, 
radial, adiabatic perturbations (pulsations), if > for 
all n. The configuration is marginally stable, if the lowest 
eigenfrequency wq = 0. 

6. Adiabatic indices 

The adiabatic index within the star, F, plays a central role 
for both the linear radial oscillations, and for the stabil- 
ity of PNS. In what follows, we will restrict ourselves to 
the case of the neutrino-opaque interior; the layer above 
the neutrinosphere contains only about 10"'^ of the to- 
tal mass, and does not influence neither the spectrum of 
radial oscillations, nor the stability of PNS. 

Under physical conditions prevailing within the hot in- 
terior of a PNS, the perturbation of a local nucleon den- 
sity, 5n, of a matter element during radial oscillations, 
takes place at constant entropy per baryon, s, and con- 
stant electron lepton number per baryon, Yj. Due to high 
density and temperature, all constituents of matter can 
be considered as being in thermodynamic equilibrium (the 
timescale of reactions leading to thermodynamic equilib- 
rium, Tioact, is much shorter, then the pulsation timescale. 



T'puis ''dyn)- The adiabatic index, governing linear per- 
turbation of the pressure within the star under these con- 
ditions, will be denoted by Fa. It will be given by 



Fa 



n 
P 



dP 
dn 



(15) 



where the derivative is to be calculated at fixed s and Yj. 
It should be stressed, that the calculation of Fa requires 
more than just the knowledge of the EOS of the PNS 
interior. Namely, one needs to know the perturbed EOS, 
under the constraint of constant s and 1;, and assuming 
the thermodynamic equilibrium of matter constituents. 

In the case of dense matter one can define several F's, 
which are physically (and numerically) different from Fa. 
The EOS of dense hot matter within a static PNS yields 
the quantity 



_n fdP 



(16) 



Strictly speaking Feos is not an "adiabatic" index unless 
the entropy s is constant throughout a star. 

In the case of a sufficiently cold matter, the reactions 
between matter constituents are so slow, that the matter 
composition remains fixed (frozen) during perturbations, 
because rj-oact ^ Tdyn- Also, neutrinos do not contribute 
then to the thermodynamic quantities of the matter. In 
such a case, the appropriate adiabatic index would read 



^ frozen — I , 

F \ an 



(17) 



where s. Ye correspond to the equilibrium model. The 
quantity Ffiozon is relevant for radial pulsations of stan- 
dard, cold neutron stars, where the difference between 
Ffrozen and Feos has interesting consequences for the sta- 
bility of the cold NS models in the vicinity of Mmax (Gour- 
goulhon et al. 1995). 

The dependence of various F's on the matter density 
in the PNS interior, for two models of the hot neutrino- 
opaque core of PNS, is displayed in Fig. 5, 6. Above nu- 
clear density (i.e., for bulk, homogeneous matter) Ffrozen 
is the highest of all F's; freezing of lepton composition 
stiffens the matter (c.f. §81 of Landau & Lifshitz (1987)). 
In the case of the isothermal core, Fig. 6, Feos is signifi- 
cantly lower than Fa and Ffrozen, for subnuclear densities 
above (this results from the fact of a significant density 
dependence of the entropy per baryon, in this density in- 
terval). For isothermal models Feos = Tstar < Ea, except 
for the region close to and below the neutrinosphere: for 
p > 10^^ g cm~^ our isothermal PNS models are convec- 
tively stable. The differences between various F's differ- 
ences are much smaller in the case of the isentropic core, 
where fixing s throughout the star and within the EOS 
makes Feos and Fa undistinguishable, except near the 
neutrinosphere, where a dramatic drop in s and in T pro- 
duces a strong deviation of Feos from both Fa and Ffrozen- 
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Fig. 6. Parameters Fa, Tbos, and Ffrozen, versus matter den- 
sity, for the isothermal hot interior with Tb = 15 MeV and with 

no trapped lopton number. The outer edge of the hot core has 
been located at = 3.5 x 10^^ gcm~^, and the rapid changes 
in the F's below results from the temperature drop in the 
neutrinosphere of the PNS. 

The differences between Ffrozen and Fa are mainly due to 
the changes of Yg throughout the star. The quantity Ffrozen 
corresponds to adiabatic pulsations keeping fixed, while 
Fa was calculated assuming beta equilibrium during adi- 
abatic oscillations at fixed Yi. It should be stressed, that 
pulsational properties of PNS arc determined essentially 
by the values of the relevant F well above the density n^. 
The outer layer with n < contains a very small fraction 
of stellar mass, and therefore rapid variations of F close 
to Til, have no effect on the global dynamics of PNS. 

7. Eigenfrequencies and instabilities 

The effects of high temperature and those of the trapped 
neutrinos influence both the EOS and the adiabatic in- 
dex of the interior of PNS. Our calculations show a 
rather strong influence of these effects on the eigenfunc- 
tions and the eigenfrequencies of radial pulsations of PNS. 
Also, these effects imply significant differences, especially 
within some (observationally interesting) interval of stellar 
masses, between pulsational properties of PNS and cold 
NS. 

In Fig. 8 we plotted the eigenfrequency of the funda- 
mental mode of the PNS pulsations, versus stellar mass, 
for several models of the EOS of the hot interior. For 



Fig. 7. Parameters Fa, Feds, and Ffrozen, versus matter den- 
sity, for the isentropic hot interior with s = 2 and with trapped 
lepton fraction Yi = 0.4. The outer edge of the hot core has 
been located at = 3.5 x 10^^ gcm""^, and the rapid changes 
in the F's below results from the temperature drop in the 
neutrinosphere of the PNS. The very steep drop in F's above 
10^** g cm~^ results from the phase transition from matter with 
nuclei to bulk (homogeneous) dense matter. 



the sake of comparison, we have presented also the cor- 
responding plot for cold neutron stars. The effect of fi- 
nite entropy and of trapped lepton number on loq depends 
rather strongly on the stellar mass. The relativistic insta- 
bility takes place very close to Mmax, so that the "clas- 
sical stability criterion" (protoneutron star models to the 
left of the maximum in the M — R plot. Fig. 4, are sec- 
ularly unstable with respect to the n = oscillations) is 
valid to a very good approximation. Generally, entropy 
and trapped lepton number soften the n = mode with 
respect to the T = models. This softening is particularly 
strong for lower stellar mass. In the case oi s = 2,Yi = 0.4 
PNS, the fundamental mode becomes secularly unstable 
at M = 0.676 Mq, very close to the minimum of the M—R 
curve. Fig. 4. The classical stability criterion gives there- 
fore a rather precise location of the instability point for the 
isentropic PNS; the models to the right of the minimum 
in the M — R curve are unstable with respect to the radial 
pulsations in the fundamental mode. However, the value of 
Mmin is dramatically larger than that for the cold (T = 0) 
neutron stars. We have M„iin[s = 2, V; = 0.4] = 0.675 Mq, 
while Mmin[T = 0] = 0.054 Mq. 
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Fig. 8. The eigenfrequencies of the fundamental mode, n — 0, 
versus stellar mass for considered models of dense matter. The 
nearly horizontal part of the curve for Tb = 15 MeV model 
corresponds to the configurations with central density below 
the value corresponding to the minimum mass of PNS. For the 
sake of comparison we also show corresponding curve for cold 
neutron stars. 

In Fig. 9 wc plotted the eigenfrequencies of the n = 
1, 2 radial modes versus stellar mass, for the isentropic 
models of the hot PNS cores. For the sake of comparison, 
we have shown also the corresponding plots for cold NS. 
The differences between the NS and the PNS n = 1, 2 
eigenfrequencies are very large, except for a narrow inter- 
val of masses in the vicity of the maximum mass. These 
differences reflect the drastic differences in the structure 
of the outer layers of the PNS and the NS models. Gen- 
erally, PNS are much softer with respect to the n = 1, 2 
radial modes. Different structure of the outer layers of stel- 
lar models with isentropic and isothermal cores, implies 
strong differences in the spectra of the lowest three modes 
of radial pulsations. Clearly, the M < 1 Mq isothermal 
models are much "softer" with respects to all considered 
modes of radial pulsations, than the isentropic ones; this 
is due to the fact that they are significantly less com- 
pact (inflated by the thermal effects) than their isentropic 
counterparts. The characteristic "turning points" and the 
intermodal intersections for the PNS models on the low- 
mass side of Fig. 9, are due to the existence of minimum 
masses. However, it should be stressed that the frequency 
intersection points, seen in the case of the case of the n = 1 
and the n = 2 eigenfrequencies, do not correspond to the 
same configuration of the PNS (while having the same 



Fig. 9. The eigenfrequencies of the first and second overtone, 
n = 1,2, versus stellar mass for considered models of dense 
matter. For the sake of comparison we also show corresponding 
curves for cold neutron stars. 



mass, these configurations have different central densities, 
see Figures 10, 11). 

The situation is less complicated in the case of the 
plots of the eigenfrequencies versus the central density of 
the stellar models, Pcentr, presented in Figures 10 and 11. 
At lower central densities the values of Wn decrease mono- 
tonically, but very slowly, with decreasing central density. 

The effects of temperature and of the differences in 
adiabatic indices F are particularly strong in the case of 
the isothermal hot interior, with T = Tao{gm)~^^'^ ■ Our 
results for the first three radial modes (n = 0, 1, 2), for 
Tb = 15 MeV, are given in Fig. 11, where we show also 
for comparison results obtained for cold NS. Three curves 
presented in Fig. 11 correspond to three different choices 
of r, governing the changes of pressure vs. density during 
pulsations. Although only results for F = Fa are physical, 
the difference between curves visualizes the role of the 
proper treatment of chemical equilibrium and adiabacity 
when the PNS star is oscillating. Feos is determined by 
the equation of state of the matter within a static model of 
PNS. As one can see from Fig. 11, the use of Feos leads 
usually to a significant underestimating of the eigenfre- 
quencies of oscillation modes. A notable exception from 
this behavior is that of the fundamental mode near the 
maximum mass of PNS, for which results obtained us- 
ing Feos are very close to those obtained using Fa- Con- 
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Fig. 10. The eigenfrequencies of three lowest radial modes, 

n = 0, 1, 2, for hot isontropic protoneutron stars with trapped 
lepton fraction Yi = 0.4, versus central density. For the sake of 
comparison we also show corresponding curves for cold neutron 
stars (thin solid lines). 




lOglo Pcentr [g Cm ^] 



Fig. 11. The eigenfrequencies of three lowest radial modes, , 
n = 0, 1, 2, for protoneutron stars with hot isothermal core 
with Tb = 15 MeV, versus central density For the sake of com- 
parison we also show corresponding curves for cold neutron 
stars (thin solid lines). 



sequently, the classical static stability criterion at Mmax 
works rather well also for the isothermal PNS. 

However, the difference between uj^'^^^ and actual loq 
increases with decreasing Pcentr and the static stability 
criterion at Mmin docs not hold. We have Afmin[7b = 
15 MeV] =0.86 Mq but configurations in the neighbour- 
hood of the mass minimum turn out to be stable with 
respect to the fundamental mode of the radial pulsations, 
calculated using the actual (physical) values of F = Fa- 

In Fig. 11 we can see an interesting phenomenon — 
abrupt, nearly stepwise changes in frequencies of the con- 
secutive (neighbouring) oscillation modes, especially pro- 
nounced for cold neutron stars, but clearly visible also 
in the case of Tb = 15 MeV, Y„ = models of PNS. 
This is the "avoided crossing" phenomenon of the radial 
modes of neutron stars. This effect is known from the 
analysis of the nonradial oscillations in "ordinary" stars 
(Aizenman, Smeyers and Weigert (1977), Christensen- 
Dalsgaard (1980)) and has been recently considered by 
Lee and Strohmayer (1996) for nonradial oscillations of 
rotating neutron stars. Our preliminary studies of this ef- 
fect indicate, that the stepwise changes of a;„ are due to 
the change of the character of the standing-wave solution 
for the eigenproblem. Namely, at the "avoided crossing' 
point the solution changes from the standing wave local- 
ized mainly in the outer layer of the star, to that localized 
predominantly in the central core. This topic will be dis- 
cussed in separate paper (Gondek et al. 1996). 

8. Discussion and conclusions 

The birth of protoneutron stars is a dynamical process, 
and a newborn protoneutron star is expected to pulsate. 
These pulsations were excited during star formation. In 
the present paper we studied linear, radial pulsations of 
protoneutron stars, under various assumptions concerning 
hot stellar core, and for a large interval of stellar masses. 

The spectrum of the lowest modes of radial pulsations 
of protoneutron stars is quite different from that of cold 
neutron stars. Generally, protoneutron stars arc signifi- 
cantly softer with respect to the radial pidsations, than 
cold neutron stars, and this difference increases for higher 
modes and for lower stellar masses. These differences stem 
from different structure of protoneutron stars, which in 
contrast to cold neutron stars have extended envelopes, 
inflated by thermal and trapped neutrinos effects. 

The standard static criteria of stability of neutron star 
models were derived under the assumption of cold, cat- 
alyzed matter (Harrison et al. 1965). We have shown, 
that to a rather good approximation, the configuration 
with maximum mass separates stable configurations from 
the secularly unstable ones (with respect to fundamental 
mode of small radial pulsations). Therefore, despite ther- 
mal and neutrino trapping effects, one can apply standard 
"maximum mass" criterion to locate the relativistic insta- 
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bility of protoneutron stars, both for the isentropic and 

the isothermal hot, neutrino opaque cores. 

The role of thermal effects increases with decreasing 
stellar mass. Static "minimum mass" criterion does ap- 
ply with a rather good precision to protoneutron stars 
with hot isentropic, neutrino opaque cores. This is due 
to the fact, that the perturbation itself conserves entropy 
(pulsations are adiabatic), and equilibration of matter is 
sufficiently fast. However, the "minimum mass" criterion 
does not apply to protoneutron stars with hot isothermal 
cores. In both cases, the minimum masses of protoneu- 
tron stars are rather large (for our models of dense hot 
matter we obtained the values of 0.675 M© and 0.86 Mq 
for the isentropic (s = 2, Yi = 0.4) and the isothermal 
(Tb = 15 MeV) models, respectively. 

Our treatment of the thermal state of the protoneu- 
tron star interior should be considered as very crude. The 
temperature profile might be affected by convection. Also, 
our method of locating the neutrinosphere was very ap- 
proximate. Clearly, the treatment of thermal effects can be 
refined, but we do not think this will significantly change 
our main results. 

Our calculations were performed for only one model of 
the nucleon component of dense hot matter. The model 
was realistic, and enabled us to treat in a unified way the 
whole interior (c;ore as well as the envelope) of the pro- 
toneutron star. However, in view of the uncertainties in the 
EOS of dense matter at supranuclear densities, one should 
of course study the whole range of theoretical possibilities, 
for a broad set - from soft to stiff - of supranuclear, high 
temperature EOS. An example of such an investigation in 
the case of the static protoneutron stars, is the study of 
Prakash et al. (1997). 

In a hypothetical scenario, proposed by Brown and 
Bethe (1994), a protoneutron star borns as a hot, lepton 
rich object, composed of nucleons and leptons. If the fi- 
nal state of cold neutron star is characterized by a large 
electron fraction, due to the appearance of a large am- 
plitude K~ - condensate (or because a large fraction of 
baryons are negatively charged hyperons), then the cold 
equation of state of neutron star matter is softer than 
that of hot, lepton rich protoneutron star. In such a case, 
A'-fbar,max(NS) < Mbar,max (PNS) , and all protoncutron 
stars with baryon mass exceeding -Mbar,max(NS) will even- 
tually collapse into black holes. Clearly, in such a case the 
problem of stability of the protoneutron star models in 
the vicinity of Mbar,max(PNS) with respect to radial pul- 
sations looses its significance. 
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